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ABSTRACT 

We report on deep spectroscopy using LRIS on Keck I and FORS2 on the VLT of a sample of 22 
candidate z ~ 6 Lyman Break galaxies (LBGs) selected by the 1775 — 2:350 > 1.3 dropout criterion. 
Redshifts could be measured for eight objects. These redshifts are all in the range z = 5.5 - 6.1, 
confirming the efficiency of the 1775 — z^q color selection technique. Six of the confirmed galaxies 
show Lya emission. Assuming that the 14 objects without redshifts are z ~ 6 LBGs, but lack 
detectable Lya emission lines, we infer that the fraction of Lya emitting LBGs with Lya equivalent 
widths greater than 20 A among z ~ 6 LBGs is «30%, similar to that found at z ~ 3. Every Lya 
emitting object in our sample is compact with r^i < Of! 14. Furthermore, all the Lya emitting objects 
in our sample are more compact than average relative to the observed size-magnitude relation of a 
large z-dropout sample (332 candidate z ~ 6 objects). We can reject the hypothesis that the Lya 
emitting population is a subset of the rest of the z ~ 6 LBG population at >97% confidence. We 
speculate the small sizes of Lya emitting LBGs are due to these objects being less massive than other 
LBGs at z ~ 6. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
starburst — early universe 



1. INTRODUCTION 

Distinguishing high redshift galaxies from interlop- 
ers at lower redshifts can be a challenging process 
(jStern fc Spinradl Il999f) . One particularly unique trait 
of high-redshift galaxies is the continuum break they 
possess as a result of absorption by neutral hydroge n 
along the line-of-sight (jMadaul 1199a iDickinsonl Il999f) . 
Putting together this feature with the particularly blue 
colors of star-forming galaxies, we have a simple but ef- 
ficient te chnique for selecting high-redshift star-forming 
galaxies (|Steidel et alJll995tlDickinsoiJll999l ). This tech- 
nique has been extensively tested through spectroscopy 
on numer ous Lyman break galaxies (LBGs) at z ~ 
2.5 - 4.5 (|Steidel et al.lll999h and an increasing number 
of galaxies and quasars at z > 4 ()Wevmann et alJll998t 
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120031 Dickinson et all 12004 IVanzella et all l2006h . In 

addition, there are now surveys that exploit narrow- 
band excess from Lya emission (see Taniguch i et al.l 
l2003h to compile la rge (> 20) samples of Lya emitters 
(LAEs) at z » 5.7 (E hQads_et_aLl |2QQ l iHu^jalJ [200l 
Shimasaku et al. 2005) and at z ss 6.6 (Kashikaw a et all 
20061: iStern et al.ll2005lh 
Several groups, utilizing these techniques, have ob- 
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ident ify ing candidate g alaxies o ut to z 
2002t iKodaira et al.l l2003t iBquwens et al 
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Taniguchi et alj 12004 IStern et all T2 005: Rh oads et all 
2004 iKneib et alJ 12004 iBouwens fc Illingworthl 120061 : 
Ive et all l2006n ~ This has led to great progress in 



our understanding of the early universe, providing 
us with con straints on both the epoch of reion- 
ization (see lLoeb fc Barkanal 120011 : IFan et all 120031 : 
iMalhotra fc Rhoads! I2004D and the ev olution of the 
global star format i on rate (SFR) density (Bouwen s et al.l 
120031 12004al 120061 : iHopkin j|2004 iGiavalisco et al.ll2004h . 
However, to ensure that these results are accurate, it 
is essential that we understand nature of the galaxy 
populations being selected by these techniques and can 
quantify their relation to the underlying population of 
star- forming galaxies at z > 6. This is paramount if we 
are to construct a self-consistent picture of the galaxy 
population from the different selection techniques. Issues 
that need to be addressed include (1) quantifying the 
distribution of Lya equivalent widths for star-forming 
galaxies at z ~ 6, (2) determining the impact of these 
results on narrowband and broadband selections, and 
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(3) using search results to construct a self-consistent 
picture of the galaxy population at z > 6. 

In this study, we utilize spectroscopy t o study the na- 
ture of z ~ 6 i775-dropouts identified by iBouwens et all 
(|2003L l2004ah in deep ACS imaging data. By contrasting 
the properties of our spectroscopic sample with similar 
selections at z ~ 3, we examine the evolution of the spec- 
tral properties of these sources with redshift. In addition, 
we pay particular attention to the implications of these 
results for future studies of high redshift galaxies using 
either narrowband or broadband selection techniques - 
noting their respective strengths and complementarity. 
A brief outline of this paper follows. A summary of the 
i-dropout sample and spectroscopic observations is pre- 
sented in §2. In §3 we present the spectra and their in- 
terpretations. In §4, we compare our results with those 
obtained with narrowband selections at z ~ 6. In §5, 
we discuss the fraction of LAEs in the i-dropout sam- 
ple, as compared with lower redshift samples. We also 
investigate a possible link between the morphologies of 
the i-dropouts and their Lya emission. Finally, in §6 we 
summarize our findings, and discuss their implications 
for future high-redshift surveys. 

In this paper we adopt the 'concordance' cosmology: 
a A-dominated, flat universe with Qm = 0-3i Oa = 0.7, 
and #0=71 km s" 1 Mpc^ 1 (jBennett et all I2003D . A ll 
magn itudes are given in the AB system (|Oke &: Gunnl 
1983). We denote the Hubble Space Telescope Advanced 
Camera for Surveys (ACS) F435W, F606W, F775W, and 
F850LP passbands as -B435, Veo6, 2775 and zsso respec- 
tively. 

2. SAMPLE + OBSERVATIONS 

2.1. Optical Broadband Dropout-Selected Survey 

Our spectroscopic sample was derived from HST ACS 
observations over two separate fields: RDCS 1252.9- 
2927 (CL1252), a cluster at z = 1 237 selected from 
the ROSAT Deep Cluster Survey (jRosati et alj 119981 
20041) . and the Ultra D eep Field parallel fields (UDF PFs; 



Bouwens et aTll2004alh A 6' by 6' mosaic of ACS Wide 



Field Camera 1775 and Z850 band images was acquired 
around CL1252, with 4' by 4' VLT ISAAC IR imaging 
of the mosaic center (|Lidman et all I2004D . The UDF 
PFs were imaged with a 4.'5 by 4f 5 ACS mosaic in the 
B435, ^606 j *775 ancl 2 850 band filters. For further details 
concerning the imaging data and our photometry over 
these fi elds, we refer the reader to IBouwens et al.1 ((20031 
l2004ah . 

The «775 and zszo band imaging data over these fields 
allo w us to readily select LBGs in the range 5.5 < z < 
6.2 (|Bouwens et al.ll2003D using an i-dropout criteria. In 
the case of CL1252, the selection was bas ed upon a ^775 — 
Z85Q > 1.3 color cut and for the UDF PFs (|Bouwens et al.l 
l2004al) . a 1775 — z$ 50 > 1.4 cut plus a null detection in Vqoq 
(2 a). Note that the colour selection we use for following 
up sources in the CL1252 field is more inclusive than th e 
*775 — z 850 > 1-5 criterion used in IBouwens et alj (2003). 
The «775 — z 850 > 1-3 criterion used for selection in the 
CL1252 field resulted from a small systematic error that 
was present in the early ZS50 ACS zeropoint. The typical 
error for the 2775 — zsso colors is 0.4 mag. In the case 
of a non-detection, the 1775 flux was set to the 2er upper 
limit. The error and upper limits are incorporated into 
all simulations that determine selection volume for the 
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Fig. 1. — An illustration of the reduction process for our two 
dimensional spectral data. The upper panel shows the unreduced 
spectrum, the middle panel shows a two dimensional fit to the back- 
ground, and the lower panel is the reduced spectrum. 



different samples. A total of 25 i-drops were identified 
in the CL1252 field, and 40 over the UDF PFs, with 10rr 
limiting magnitudes of Zs,bo = 27.3 and 27.8 respectively. 
This corresponds to surface densities of 0.5 ±0.1 objects 
arcmin -2 for CL1252, and 1.4 ± 0.2 objects arcmin -2 
for the UDF PFs. 

In Table [T] and Table [H the objects satisfying the i- 
dropout criteria for spectroscopic follow-up are listed, 
along with their zsso-band, J-band and i^s-band magni- 
tudes, and 1775 — Z850 colors. A total of 22 objects were 
observed spectroscopically (14 in CL1252 and eight in the 
UDF PFs). These objects were randomly selected from 
the entire j-dropout sample for spectroscopic follow-up, 
except for the most luminous dropout in the CL1252 field 
(BD38) which was explicitly targetted for follow-up. In 
addition, i-dropout galaxies were the primary targets for 
both the CL1252 and UDF PFs spectroscopic observa- 
tions. The faintest galaxy on our mask had a zsso-band 
magnitude of 27.5. 

2.2. Keck I LRIS 

We us ed the Low Re solution Imaging Spectrograph 
(LRIS) (|Oke et al.lll995l) on the Keck I 10m telescope 
for spectroscopic follow-up of all 22 objects in our spec- 
troscopic sample. A total integration time of 16200 s for 
CL1252 objects and 7200 s for UDF PFs objects was ob- 
tained using the 600 line mm -1 grating blazed at 8500 A 
with a slitlet width of 1" and a minimum slitlet length of 
10". The scale (AA res ) of this configuration is 1.28 A per 
pixel. A series of eight dithered exposures was taken to 
aid in cosmic ray removal, and to enhance removal of the 
fringing in the near-IR region of the spectra. The tele- 
scope was offset l"-3"along the slits between exposures, 
which ranged from 1200 to 2250 s in duration. The ob- 
servations were taken under photometric conditions on 
the nights of UT 2004 February 13-14. 

All data reductions were conducted using a slit mask 
reduction task developed by Daniel Kelson. This 
task yields cleaner background subtraction and cosmic 
ray removal t han standard slit spectroscopy procedures 
(|Kelsonll2003D . By performing the sky subtraction before 
the data have been rectified, this methodology minimizes 
reduction errors. One example of this reduction proce- 
dure is presented in Figure[T] Wavelength calibration was 
performed using the night skylines. The spectra were flux 
calibrated using a sensitivity function derived from obser- 
vations taken t he same night of the spectrophotometric 
stand ard HZ44 (|Massev et alj|1988t iMassev fc Gronwalll 
1990). 



Spectroscopic Confirmation of z ~ 6 LBGs 
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TABLE 1 

Properties of our RDCS 1252-2927 Spectroscopic Sample 
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Previous ID a 




RA 






Dec 
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BD58 


1252- 


2585-3351 


12 


52 


52 


283 


-29 


28 


01 


74 


25 


7±0 


1 


25.6±0 


3 


25.4±0.3 


2.0 





20 






BD44 


1252- 


5058-5920 


12 


53 


01 


745 


-29 


26 


03 


87 


25 


9±0 


2 


26.3±0 


4 


25.5±0.4 


>1.8 





19 






BD46 b - d 






12 


53 


05 


424 


-29 


24 


26 


22 


26 


1±0 


1 


NA 




NA 


1.6 





11 


5 


914 


BD57 






12 


52 


49 


441 


-29 


27 


53 


24 


26 


1±0 


2 


25.7±0 


3 


24.9±0.2 


1.5 





19 






BD00 d 






12 


52 


42 


927 


-29 


29 


20 


05 


26 


1±0 


1 


NA 




NA 


1.3 





11 


5 


942 


BD66 






12 


52 


54 


624 


-29 


24 


56 


18 


26 


1±0 


1 


NA 




NA 


1.4 





13 






BD40 






12 


52 


59 


511 


-29 


26 


58 


13 


26 


3±0 


1 


26.2±0 


3 


25.9±0.4 


1.3 





15 






BD22 






12 


52 


53 


346 


-29 


27 


10 


58 


26 


5±0 


2 


NA 




NA 


1.4 





13 






BD62 


1252- 


3729-4565 


12 


52 


56 


746 


-29 


27 


07 


63 


26 


7±0 


2 


26.3±0 


3 


>26.2 


>1.8 
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BD48 


1252 


-3497-809 


12 


52 


42 


754 


-29 


27 


18 


89 


27 


0±0 


2 


NA 




NA 


>1.6 





11 






BD36 






12 


52 


54 


944 


-29 


25 


57 


52 


27 


1±0 


1 


26.3±0 


4 


25.2±0.4 


1.5 
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Note. — Right ascension (hours, minutes, seconds) and declination (degrees, arcminutes, arcseconds) use the J2000 equinox. All 
magn itudes given are AB. An entry of NA indicates that the object was outside the ISAAC coverage we had on CL1252 dLidman et all 
2004). Each of these objects were observed for four hours with LRIS. 

a IBouwens ct al. (2003). ^ FORS2 22.3 hr integrated spectra also obtained. Lacking detectable Lya emission; absorption line 



redshift. Lya emission detected; emission line redshift. 



TABLE 2 

Properties of our UDF Parallel Spectroscopic Sample 



Object ID 


Previous ID a 


RA 


Dec 


2850 


«775 — 2850 




z 


UDF PFsl iO 




03 32 35.604 


-27 57 37.51 


25.5±0.1 


1.9 


0.25 




UDF PFsl il 




03 32 34.892 


-27 57 10.74 


25.7±0.1 


1.5 


0.11 




GOODS i6 b 




03 32 39.803 


-27 52 57.88 


26.1±0.1 


1.4 


0.10 


5.540 


UDF PFs i4 b 


UDFP1-3851-2438 


03 32 43.959 


-27 56 43.87 


26.3±0.1 


1.6 


0.11 


5.857 


UDF PFs i9 


UDFP1-4650-3354 


03 32 40.707 


-27 57 24.28 


26.4±0.1 


1.7 


0.16 




UDF PFs il b 


UDFP1-2954-1152 


03 32 48.368 


-27 55 54.82 


26.9±0.2 


1.5 


0.09 


6.005 


UDF PFs i2 b 


UDFP1-2309-1628 


03 32 46.410 


-27 55 24.32 


27.1±0.1 


>2.5 


0.10 


6.083 


UDF PFs i0 c 


UDFP1-3407-1028 


03 32 48.952 


-27 56 16.97 


27.5±0.2 


1.7 


0.10 





from an alignment star spectrum. Therefore it is excluded from the sample. 



2.3. VLT FORS2 

Three CLf252 i-drop objects (BD38, BD03, and 
BD46) were also observed with the Focal Reducer/low 
dispersion Spectrograph 2 (FORS2) on the 8.2-m VLT 
YEPUN Unit Telescope in Mask Exchange Unit mode. 
These objects were observed as part of a large observing 
program aimed at acquiring very deep spectra of both 
field and cluster elliptical galaxies in the redshift range 
0.6-1.3. We used the 600z grism with the OG590 order 
separation filter, which yielded a scale element of 1.64 A 
per pixel. The data presented here were taken from UT 
2003 February 28 through March 2 and April 27 through 
May 24 with a median seeing of 0'.'65 and clear condi- 
tions. The observations were carried out in a series of 
four dithered exposures with equal exposure times rang- 
ing from 14 to 30 minutes each, yielding a total integra- 
tion time of 22.3 hours for the sources. All exposures are 
added and weighted such that optimal signal-to-noise ra- 
tios were obtain ed. Details concern i ng th e data reduction 
are provided in Ivan der Wei et ail ((20051) . 

3. RESULTS 

3.1. Keck Data 

Six emission line objects were detected out of the ob- 
served 22 object sample. No other spectral features were 
found (including continuum) in the remaining spectra. 
Figures [2^ - f show our spectra of these objects, with 



lines evident at 8440 A, 8404 A, 8514 A, 8609 A, 8335 A, 
and 7950 A. Since O n A3727 would be resolved into a 
doublet structure with our spectral resolution, this inter- 
pretation for these lines could clearly be rejected. More- 
over, line interpretations as H/3 A4861.3, Ha A6562.8, 
and O in A5006.8, 4958.9 could also be discarded due to 
the lack of nearby lines and the 1775 — z^q >1.3 colors 
(strong continuum breaks) of our objects. For objects 
BDOO, BD46, UDF PFs il, and UDF PFs i2, the asym- 
metry of the line profiles, with absorption on the blue 
side, is consistent with absorption by a thick Lya forest. 
Therefore, it seems quite clear that these emission lines 
are Lya A1215.67, placing these objects at z — 5.942, 
5.9214, 6.005, and 6.083, respectively. For UDF PFs i4 
and GOODS i6 0, the emission features are heavily con- 
taminated by sky lines, making it difficult to ascertain 
their level of asymmetry and thus whether they are likely 
Lya. We tentatively identify these emission features as 
Lya placing these objects at z — 5.857 and 5.540, re- 
spectively. However, in the discussion which follows, we 
shall also consider the implications that these objects are 
null-detections or low-redshift contaminants. 

To quantify the asymmetry of the four uncontaminated 
emission lines, we adopted the two asymmetry parame- 
ters developed by Rhoads et al. 2003. The "wavelength 
ratio", a a, and "flux ratio", a/ are defined as, 
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Alfj, r ~ A p 
An A 



10.6 



a f 



(i) 



(2) 



respectively. Here, A p is the wavelength of the peak 
of the emission, and Aio,& and Aio, r are the wavelengths 
where the flux exceeds 10% of its peak value on the blue 
and red side of the emission. The resulting values, shown 
in Tabled range from a w = 1.8-2.6 and a/=1.4-2.4. This 
is in agreement with va lues determined fro m Lya emit- 
ting objec^^z^^^^^s^MlMM^^j and at z ~ 
5.7 (jRhoads et alJ [2003). These asymmetry parameters 
are inconsistent with those of the z ~ 1 O n A3 727 dou- 
blet (a w < 1 and a f < 1; iDawson et all (|2004h ). Note 
that for object BD46, the FORS2 spectrum (see §3.2.1) 
was used to calculate asymmetry parameters. 

The spectrum of UDF PFs iO was contaminated by 
bleeding from an alignment star, impeding detection of 
Lya emission or a continuum. Thus, this object is ex- 
cluded from our sample in the analysis below. 

The two objects UDF PFs il (z 850 = 26.9) and UDF 
PFs i2 (zs5a = 27.1) arc among the faintest z > 5 ob- 
jects in the zsso band to be selected by the Lyman break 
technique with confirmed redshifts, very similar to the 
objects GLARE 3001 and 3011, with z 850 = 26.37 ± 0.0 6 
and 27.15 ± 0.12 respectively (|Stanwav et alJ l2004bT >. 
Furthermore, UDF PFs i2 is faintest object in our sam- 
ple (excluding UDF PFs iO), and is the object with the 
highest confirmed redshift. 

Table [3] gives the measured redshifts, Lya fluxes 
(/iya), rest-frame equivalent widths (WLy a )> and 
FWHMs of the six Lya emitting LBGs. We measured 
the flux in the emission lines between the red dotted lines 
illustrated in Figures 2 and 3. In some cases, there were 
residuals from the sky subtraction that could contami- 
nate the line flux. We estimated this additional contribu- 
tion from the sky lines by measuring the residual flux in 
two-dimensional sky-subtracted spectrum outside of the 
extraction region used for the one-dimensional spectrum. 
All values quoted were measured directly from the spec- 
tra, and do not account for flux lost due to HI gas absorp- 
tion in the objects. Therefore, the quoted fluxes are lower 
limits. The observed emission fluxes from our Lya emit- 
ting sample are very similar to those found for other z ~ 6 
objec t s selected by the i-dro pout method (|Bunker et al.l 
120031 IStanwav et all l2004afbl ). and for z « 5.7 and z 
6.6 Ly a emitters selected using narrowband filters 
dHu et alJ 12004 [Rhoads et alj|2003t iLehnert fc Bremerj 
2003tlAiiki eta l. 2003; Kurk et al. 2004; Tani guchi etafl 
2004 t IStern et al.ll2005l : iNagao et alJl2004fc iRhoads et all 
2004J). To determine the continuum for each spectrum, 



we assumed it satisfied a power law of the form 
/(A) = / o (A/0.9 M m)^ 



(3) 



with (3 = —1.1, which is the average of 198 Lya emitting 
LBGs at z ~ 3 dShaplev et alJ l2003L f„ was determined 
by fitting the zsso magnitude. The equivalent widths 
of the Lya emissions were calculated directly from these 
estimated continua. These widths were not very sensitive 
to the assumed value of (3. 



3.2. VLT Data 

3.2.1. BD46 

The two-dimensional FORS2 spectrum of BD46 
showed a robust emission line at the central wavelength 
of 8406 A. Two- and one-dimensional extracted spectra 
for this emission are shown in Figure [3^. A clear asym- 
metric profile for the line is seen, which we identify as 
Lya yielding a redshift of z =5.9 14 for this source. Fur- 
thermore, a weak continuum was detected redward of 
the emission, with / = 3.4 ±1.0 counts A -1 , shown in 
Fig. Wp- One item to note is BD46's Lya emission was 
detected with LRIS before the VLT data was available. 
Our methodology accurately measured this redshift (it 
was within Az = 0.0005 of the FORS value), despite a 
much weaker signal due to the smaller integration time 
in our Keck LRIS data. 

3.2.2. BD03 

While this object was undetected in the LRIS (4 hour 
integration) spectrum, the FORS2 spectrum shown in 
Figure H contains a noisy but clearly flat continuum red- 
ward of 8500 A, with a sharp discontinuity at 8485 ± 
3 A, where the flux drops to approximately zero blue- 
ward of 8400 A. Using a constant step function to fit 
the continuum, we find the average flux density above 
the discontinuity is / red (8500 - 9100 A) = 5 ± 1 counts 
A" 1 . Below the break, this reduces to / Wue (7700 - 8300 
A) = 1 ± 1 counts A -1 , consistent with no detected flux. 

This discontinuity is much larger than would be ex- 
pected from UV/optical spectral breaks of galaxies at 
rest wavelengths 4000 A [D(4000)] 2900 A [B(2900)] 
and 2640 A [B(2640)] ([Stern et alJ l2000t ISpinrad et all 
Il997l[l998h . Another option is that this object is a radio- 
loud broad absorption line quasar. Such objects can have 
large continuum breaks of this amplitude. Ho wever, such 
objec ts also possess very red near-IR colors (Ha ll et all 
11997f t. unlike BD03 (J - Ks < -0.03). Furthermore, the 
resolved morphology of BD03 is atypical for luminous 
AGNs at high redshift, which are generally unresolved. 

The abruptness of the break and the lack of slope in 
the continuum at larger wavelengths are strong signs that 
the discontinuity is due to absorption by intergalactic HI 
gas. Coupling this with the 2775 — Zsm > 2.1 colors (the 
second reddest in our spectroscopic sample), we interpret 
this object as being a starburst galaxy at z = 5.98 ± 0.1. 
Using the LRIS spectrum we place an upper flux limit of 
8.0 x 10~ 19 ergs cm -2 s _1 on the Lya emission. 10 

3.2.3. BD38 

This object was also undetected in the LRIS obser- 
vations, but yields a strong continuum in the FORS2 
spectrum. Figure [5] shows the two- and one- dimensional 
spectra of the object, along with the ACS and ISAAC 
imaging. A precipitous continuum break is clearly seen 
at w 7900 A, reducing the continuum from 27 ± 1 to 
± 1 counts A -1 . Also, several strong absorption features 
are found at 8205.77 A, 8485.38 A, 8498.43 A, 8699.72 



10 All flux limits are for a 3 a detection e xtracted over 4 A 
(« 3AA res ). Note that IStanwav et al.l (|2004afl quoted a 10.5 hr 
DEIMOS spectra flux limit of 2.0 X IO -18 ergs cm -2 s -1 , which 
is for a 5 a detection extracted over 8 A. 
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TABLE 3 

Objects With Detectable Lyq Emission 



Object ID 


z 


f Lya ' 


W Lya h ' c 


FWHM C 


a w /a,f d 






(10 — 17 ergs cm -2 s _1 ) 


(A) 


(A) 


BD46 


5.914 


0.90 


24 


6 


1.9/1.8 


BDOO 


5.942 


3.50 


150 


10 


1.8/1.4 


GOODS i6 


5.540 


1.34 


31 


19 


UDF PFs i4 


5.857 


0.70 


34 


10 




UDF PFs il 


6.005 


1.10 


65 


13 


2.6/2.4 


UDF PFs \1 


6.083 


0.96 


64 


10 


2.3/1.5 



a Flux from Lya emission. Measured rest-frame equivalent width. c Not corrected for 
IGM absorption.*^ Emission line asymmetry parameters ( Rhoads ct al. 2003) 



A, 9079. 21 A, and 9f40.84 A in the spectrum. Associ- 
ating the break with the UV optically thick Lya forest, 
and the absorption lines with the strong interstellar ab- 
sorption features typical of LBGs at l ower redshifts (see 
IShaplev et ail 120031 : lAndo et al.ll2004h . yields a redshift 
determination of 5.515 ± 0.003. 

BD38 is an unusually bright i-dropout. Its zsso-band 
magnitude is 24.3, making this object more than 1.3 mag 
brighter than any other i-dropout in our spectroscopic 
or photometric sample. In fact, to our knowledge, this 
object is still the brightest z ~ 6 LBG discovered to date 
with a confirmed redshift. The complete spectrum, along 
wi th a more thorough a n alysis of this object is presented 
in lDow-Hvgelund et alJ (|2005[ ). 



4. DISCUSSION 
4.1. Confirmation and Contamination 

Of the total 21 galaxies in our spectroscopic sample 
(excluding UDF PFs iO), we identified spectral features 
for eight objects. Of these eight, we confirmed that six 
i-dropouts from our sample are z ~ 6 galaxies, with the 
strong possibility that two more are. Therefore, the suc- 
cess rate of our survey is assuredly 29%, and most likely 
38%. For the smaller integration times obtained with 
LRIS (4 hours), the completeness drops to 25% (19% 
disregarding GOODS i6 and UDF PFs 14). The confir- 
mation rate in our longer 22.3 hr FORS2 integrations is 
100% (three out of three). Interestingly, the two reddest 
*775 — ^850 color objects, the two brightest zgso objects, 
and the faintest zsso object in our sample are confirmed 
z ~ 6 galaxies. 

There was no evidence from the spectra that any of 
the objects in our sample are low-redshift contaminants 
despite the long integration times. This is consistent 
with these objects being z ~ 6 continuum objects with- 
out strong Lya emission. However, if we associate the 
emission lines of GOODS i6 and UDF PFs i4 with 
those of low redshift interlopers (for which we have no 
evidence) the contamination of the sample could be as 
high as 14%, though this is unlikely. Support for this 
point is pr ovided by the low resol ution spectroscopy con- 
ducted by iMalhotra et all (|2005l) using the ACS grism 
to test the i-dropout selection. Using a liberal color cut 
°f *775 ~ z 850 > 0.9 (versus our 1775 — 2:350 > 1.3 crite- 
rion), they verified that 23 out of 29 candidate i-dropouts 
showed a strong continuum break indicative of galaxies 
at z ~ 6. In addition, they found that a colour cut 
of ?775 — ^850 > 1.3, while suffering from modest (20- 
30%) incompleteness, is subject to a contamination rate 



of only 7% (one out of 14). 11 Since the estimated con- 
ta mination in the l arge (506 object) i-dropout selection 
of iBouwens et all (2006) is similarly small (<8%), the 
present findings are perhaps not too surprising. This 
lends strong support to the notion that the objects in 
our sample without detectable features are primarily un- 
confirmed z ~ 6 galaxies, and are not low redshift inter- 
lopers. 

To see if any of the unconfirmed galaxies in our CL1252 
spectroscopic sample are low-2 cluster members, we 
coadded the LRIS spectra associated with these sources, 
in an effort to increase the signal to noise of the obser- 
vations. If these objects were predominately z — 1.23 
early-type contaminants, a break at «8800 A would ap- 
pear. No continuum break was detected in the coadded 
spectra, nor was any excess continuum found. Therefore, 
we concluded that the majority of the objects in our sam- 
ple are not cluster members. We should note however 
that this null detection is not very surprising considering 
that BD38 and BD03, which are the brightest objects 
in our sample and are confirmed z ~ 6 objects via the 
FORS2 observations, were undetected in the LRIS ob- 
servations. The main conclusion of this analysis is that 
deeper spectroscopy is necessary to ascertain the nature 
of these galaxies. 



4.2. Star Formation Rate 

Aside from Lya A1215.67, no other emission lines 
were found in our sample of galaxies with Lya emis- 
sion. Another sometimes strong line in our wavelength 
range is the rest-UV line N v A1240. This emis- 
sion feature is strong in active galactic nuclei (AGNs), 
with line ratios f rom a composite quasar spectra of 
(hy a /fav = 4.0) (|Osterbrocklll989h . The flux limits at 



rest-frame 1240 A are / NV < (0.8, 0.8, 1.0, 2.4, 2.2, 1.0) x 
10 -18 ergs cm -2 s _1 , implying lower limits of the 
fhya/hv > (11.3,44.8,13.4,2.9,5.0,9.6) line ratios. 
Furthermore, no object in our sample was detected 
in Chandra or XMM-Newton X-ray imaging data 
(|Rosati et all 12004 iGiacconi et alJ l2002f >. The lack of 
any N v A1240 emission or X-ray detection suggests that 
the Lya emission is not produced by an AGN. Thus, we 
believe the Lya photons are primarily from hot young 

11 This contaminant was a star. The methodology of 
IBouwens eFaLI (f2003. !2004aL [20061) uses the SExtractor stellarity 
parameter IjBertin k A rnouts 1996) to exclude stars from our i- 
dropout selection. Hence, this contaminant would likely have been 
rejected from our spectroscopic sample. Thus, o ur null contami- 
nation results are consistent with the findings of Malhotra ct al. 
p005h . 



() 



stellar populations consistent with the Lya emitting ob- 
jects being starbursting galaxies. 

Now we estimate the star formation rates of our eight 
confirmed and 13 unconfirmed objects in our z ~ 6 sam- 
ple. Two different methods are used, the first relying the 
Lya emission flux (SFRr, ya ) to make this estimate, and 
the second relying on the [/U-continuum flux (SFRuv)- 
The relationship between the Lya flux and the SFR of a 
galaxy is given by 



SFR 



Lya 



9.1 x 1(T 43 L Lya M Q yr" 



(4) 



where Lu ya is the Lya luminosity in units of ergs 
s _1 , assuming the Salpeter initial mass function with 

{^Tljoweri Slower 

) = (0.1 Mq, 100 M Q ) (|Salpeterlll955l : 
iBrocklehurstl I197H iKennicuttl IT998) . The SFRs deter- 
mined by this method yield lower limits, due to absorp- 
tion of Lya photon s by dust grains within the galaxy and 
by the Lya forest (|Hu et al.ll200l ). The SFRs can also 
be derived from the UV continuum luminosities (Luv) 
at A = 1500 A, and using the following relation 

f<S> yr" 1 (5) 



SFRuv = 1-4 x 10~ 28 L„v M, 



wher e Luv is in units of ergs s 1 Hz 1 (jMadau et alJ 
1998). For the nine objects with ISAAC imaging, the 
slope of the continuum was derived from the zsso, J ■, 
and Ks magnitudes assuming the UV spectrum can be 
described by f\ oc A' 3 . The continuum slope used for 
the remaining objects was the consistent value of /3 = 
-2.0 found from ffg T/NICMOS imaging of 26 i-dropout 
objec ts in the UDF (jStanwav et al.ll2005t iBouwens et al.1 
2006). We adopted a different /3 than that used in §3.1, 
because of the obs erved /3 dependence on Lya emission 
strength at z ~ 3 (iShaplev et al.|[2003h . For the objects 
in our sample for which we cannot measure redshifts, 
the mean redshift for z-dropouts in our selection, z = 5.9 
(|Bouwens et al.ll2~003l ). is assumed. 

For objects in the C L1252 field, we used the results of 
lLombardi et al.l (|2005l ) to correct for gravitational ampli- 
fication by the cluster potential. For each object, we as- 
sumed a Non-Singular Isothermal Sphere lensing model 
with a free central position. This yields a best-fitting 
velocity dispersion of 1185 km s . Again, for the non- 
detected objects the mean redshift z = 5.9 was used. 
The resulting gravitational magnifications are presented 
in the second column of Table |U Most objects are two 
or more Einstein radii away, being only weakly lensed. 
However BD22 and BD62 are close to the Einstein ra- 
dius, and therefore the magnifications quoted for these 
sources are highly uncertain. Hereafter only delensed 
values for objects in the CL1252 field are quoted, except 
for BD22 and BD62 where observed values are retained. 

In Table 2] we present Lu ya , Luv, SFRr, ya , and 
SFRuv for the Lya emitting objects in our sam- 
ple. Errors in Luv and SFRuv are derived from the 
zsBo-band flux errors. As is seen elsewhe re (jHu et al.1 
12001 iKodaira et all 120031: lAiiki et al.ll2003h . for the five 
Lya emitting objects the ratio SFRu ya / SFRuv ranges 
widely from 27% to 127%. This is m ost likely due t o 
absorption of Lya photons by the IGM (|Hu et al.| [2002). 

4.3. Selection Completeness for Lya Emitters 

Lya A1215.67 emission is by far the most prominent 
spectral feature we can use for redshift determination of 
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Fig. 6. — Observed (upper) and delensed (lower) 1775 — 2850 
colors versus zgso-band magnitudes for our spectroscopic sample. 
The black squares are sources where no spectral features were found 
red circles are the sources with Lya emission, and blue triangles 
are sources with no detectable Lya emission but for which we can 
derive a redshift (BD03 and BD38). Objects undetected in the 1775 
band are shown at their 2cr lower limits and include upward arrows. 
The sources in CL1252 that are not delensed are shown with right 
arrows. The upper left orange error bars are typical for objects 
detected in the 1775 filter. Our 1775 — 2850 > 1-3 selection criterion 
is denoted with the blue dashed line. Aside from UDF PFs i2, 
the distribution of Lya emitters is constrained to zsso > 26.1 and 
i-775 — z 850 < 1-6 (demarcated by dotted lines). In the observed 
plot almost half ( 9 out of 21 ) lie outside this region. However, this 
trend is absent after delensing the objects. Lya emitting objects 
likely have bluer 2775 — zsso colors than other i- dropouts as a result 
of the contribution of Lya flux to the i77$-band magnitudes (see 
Figure 7). 



z > 5 objects. When the Lya emission is strong, it can 
have an effect on the observed broadband fluxes. At z 
< 6, Lya flux falls in the i-band, yielding i — z colors 
that are smaller than without such emission. This effect 
is notable in Figure [6l where we show the ^775 — zsso 
colors vs. zsso-band magnitudes of our sample. Aside 
from UDF PFs i2, which is at z = 6.083, the Lya 
emitting objects inhabit the re gi on ^7 75 — z 850 < 1.6. 
Moreo ve r, from iVanzella et al.l (12006D: iVanzella et al.l 
(120051); iDickinson et al l (l2004fk iNagao et al.l (I2004D: 
IStanwav et alJ (|2004aU H): iBunker et all (|2003l) and this 
work, only two confirmed Lya emitting objects have 
been detected by the i-dropout technique for z < 5.78. 
Once Lya redshifts into the Zs^q band, the decrement 
can be larger than it would be without emission, such 
as for the z = 6.33 Lya emittin g z-dropout object S DF 
J132440.6+273607 discovered bv lNagao etaTl (|2004h . 

In Figure [7j we present the redshift distribution of 
spectroscopically confirmed z ~ 6 i-dropouts. The 
uppe r histogram c ontains o bjects from th is wor k 
and IStanwav et "all (|2004allbh : iMalhotra et~aT1 (|2005h : 
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TABLE 4 

Star Formation Rates for our z ~ 6 Sample 



Object ID /x a LLj, Q b iyv c SFR Lya d SFR uv e 

(10 42 h^ 2 . ergs s" 1 ) (10 28 ergs s" 1 Hz" 1 ) (h^ M Q yr" 1 ) (/^ 2 M yr" 1 ) 



BD38 


1.3 


<0.5* 


27.1±2.7 


<0.5 


38.1±3.8 


BD03 


1.2 


<0.3 f 


lO.Oil.O 


<0.4 


13.4±1.3 


BDOO 


1.1 


13.24 


6.8±0.7 


12.1 


9.5±1.0 


BD46 


1.1 


3.37 


6.7±0.7 


3.1 


9.5±1.0 


GOODS i6 




4.31 


6.7±0.7 


3.9 


9.4±1.0 


UDF PFs i4 




2.56 


6.1±0.6 


2.3 


8.5±1.0 


UDF PFs il 




4.27 


3.5±0.7 


3.9 


5.1±1.0 


UDF PFs i2 




3.84 


3.1±0.3 


3.5 


4.3±0.4 


UDF PFs iO 




NS 


12.9±1.3 




18.0±1.8 


UDF PFs-IDROPl 




NS 


12.9±1.3 




18.0±1.8 


BD27 


1.5 


NS 


7.6±0.8 




10.6±1.1 


BD44 


1.2 


NS 


7.3±1.5 




10.1±2.0 


BD58 


1.6 


NS 


6.6±0.7 




9.3±1.0 


BD66 


1.2 


NS 


6.2±0.6 




8.67±0.9 


UDF PFsl i9 




NS 


5.6±0.6 




7.9±0.8 


BD57 


1.3 


NS 


5.3±1.1 




7.7±1.5 


BD22 


4.22 


NS 


5.1±1.0« 




7.2±1.4« 


BD62 


35.1 


NS 


4.1±0.8 g 




6.2±1.2« 


BD40 


1.6 


NS 


3.9±0.4 




5.4±0.5 


BD48 


1.2 


NS 


2.7±0.5 




3.8±0.8 


BD36 


1.4 


NS 


2.4±0.2 




3.4±0.3 



Note. — Uncertainties in Luv and SFB.uv ar c derived from the Zgso-band flux errors. All quoted luminosities and star 
formation rates for CL1252 cluster objects arc corrected for possible lensing, except for BD22 and BD62 where the observed 
values are used. An entry of NS indicates that no spectral features were detected for the object. Where no redshift could be 
determined, the quot ed UV luminosities arc for the object at the mean redshift z — 5.9 estimated for z-dropout selections 
IBouwens et al.ll2003l) . 

a Gravitational magnification due to the CL1252 cluster potential^ Lya luminosity. UV continuum luminosity at A=1500 
A.^ SFR derived from LL yo .. e SFR derived from Luv? No detectable Lya emission. s Not corrected for gravitational lensing 
and thus highly uncertain. 

IVanzella et alJ (|2006l ). The lower histogram only con- 
tains those objects which have measurable Lya emission. 
In addition, the ex pected distribut i on of z'-drop outs from 
the simulations of IBouwens et al.l ((2003, 200q) are plot- 
ted. From these curves it is evident that as the strength 
of Lya emission increases, the mean redshift of the i- 
dropout population increases quite noticeably, while the 
width of this distribution narrows somewhat. Hence, 
there is a substantial bias against lower redshift (z < 
5.7), strong {Wi, ya > 50A) Lya emitting i-dropouts be- 
ing selected as i-dropouts. However, this effect is small 
for weaker {Wl V u =25 A) Lya emitting objects. Consid- 
ering that over half of the Lya emitting objects in our 
sample have W^ya <30 A, this should only result in a 
modest bias. In § !5.1l we discuss this issue further. 

To understand how Lya emission affects z-dropout se- 
lection and how this influence evolves with redshift, we 
modeled the spectra of the Hu et al. (2004) (hereafter 
Hu04) sample of 19 confirmed z « 5.7 LAEs, all identi- 
fied by a narrowband selection. We plot the colors for 
Hu04 in Figure[H Spectra of these objects were assumed 
to satisfy the power law given in Eq. [3] The flux decre- 
ment due to the IGM was modeled as a simple trans- 
mission coefficient, where initial values fo r the fit were 
obtained from the tables of Songaila (2004|). The contin- 
uum level was set to zero blueward of the Lyman limit 
A912. The IGM transmission, f a , and total Lya emission 
integrated flux were varied as three fit parameters until 
the quoted NBgi^Q, Cousins Z , and Cousins / magnitudes 
in Hu04 were reproduced using the IRAF software pack- 
age Synphot. These fiducial spectra were then artificially 
rcdshifted or blueshifted across the space probed by the 
«775-dropout technique in intervals of Az = 0.1. At each 
interval, the 1775 and £850 magnitudes were recomputed 



using the model for each galaxy. 

By directly applying our search criterion on these cal- 
culated 1775 and Z850 magnitudes, we can determine the 
efficiency of the «775-drop method at selecting LAEs at 
various redshifts. Using the relation 




7 

= ^p e (zi)AV(^) (6) 

i=l 

where n is the expected value of LAEs, p(z) is the vol- 
ume density of Hu04 LAEs at redshift z, p is the true 
volume density of Hu04's objects, e(zi) is the fraction of 
Hu04 LAEs we would select at redshift z\ (i.e., selection 
probability), AV(zi) is the volume element for the red- 
shift interval z\+\ — Zj, and the sum ranges from z\ = 5.5 
to Z7 = 6.1. This relation yields an expectation value of 
3.0 strong LAEs in the CL1252 field. 

This expectation value is an upper limit due to two 
additional effects. First, the probability of a z ~ 6 ob- 
ject being selected by the 1775 — zsso dropout technique 
is a strong function of reds hift, dropping off rapidly a t 
the highest redshifts probed (|Bouwens et al.|[2003ll2006l) . 
This is due to surface brightness selection effects at work 
at t he high-redshift end o f our 2^6 i-dropout selection 
(see IBouwens et aT1l2006l ). Second, Lya emission red- 
shifted into the OH bands will be much more difficult to 
detect than emission in the region between the bands, 
and so the comoving volume probed will be smaller than 
we assumed. Mitigating this fact is the quoted flux limit 
of Hu04 is 2 x 10~ 17 ergs cm -2 s _1 . This is slightly 
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7. — Observed distribution of spectroscopically confirmed 
i-dropouts (2775 — 2850 > 1-3J. The upper histogram 



(2004b ), 



this work and 


Stanwav et al.\ 


Malhotra et al. 


S200SI), and 



I200id). \Stanwav et al.\ 

I Vanzella et al. (2006). The lower histogram contains objects with 
detectable Lya emission. Overplotted are simulated redshift dis- 
tributions for four different Lya equivalent widths Whya, indicat- 
ing the impact of these equivalent widths on the redshift distribu- 
tion of the i-dropout selection. Overall, the simulations respectably 
encompass the observed distribution, and suggest that most star- 
forming galaxies at z ~ 6 only have modest Lya equivalent widths. 



greater than the average flux limit for our CL1252 spec- 
tra for a Lya emission in the redshift interval probed (1.9 
x 10~ 17 ergs cm~ 2 s _1 ). 

Hence, we believe the expectation value of 3.0 LAEs 
to be consistent with our spectroscopic results of two 
Lya emitting LBGs (BDOO and BD46) in our CL1252 
sample of 12 objects. Moreover, it appears our survey 
is complete in identifying objects that have measurable 
Lya emission. Coupling this with our likely minimal low- 
rcdshift contamination rates (see §4.1|) suggests that Lya 
emitting LBGs only represent a fraction of i-dropout se- 
lected z ~ 6 LBGs. 

We estimated the total number of expected LAEs over 
the redshift range of our i-drop out selection as s uming 
the Hu04 volume density and the iBouwens et alj (|2003l ) 
covolume. This yields a strict upper limit of 5.6 expected 
LAEs in our CL1252 sample. From this estimate, the i- 
dropout selection misses at most 2.6/5.6, or ~46% of the 
total LAE population, though most of those missed are 
at the low redshift end of the i-dropout range. 

5. TRENDS IN Z ~ 6 LYMAN BREAK GALAXIES 

5.1. Fraction of Lya Emitting Galaxies 

In our sample, we find six objects (four if GOODS i6 
and UDF PFs i4 are disregarded) with rest-frame equiva- 
lent widths strong enough to be detected as narrowband 
excess objects (W Lya >20 A). Thus 29% (19%) of our 
sample (excluding UDF PFs iO) would be recovered by 
narrowband surveys. Similar results of 33%, 33% and 
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Fig. 8. — i-z vs. NBgiso-z color distribution of two narrowband 
selected samples. The red circles are spectroscopically confirmed z 
~ 5. 7 Lya emitters and the black filled squares are unconfirmed 
objects from \Hu et al.\ 1200 A) . The green open squares are na rrow- 
band excess candidate Lya emitters from VAiiki et alj If MOOS) . No 
excess flux in the zgso-band was measured for two objects in this 
sample; hence they are excluded from the plot. A typical i-dropout 
color selection of 1775 — zg^o > L3 is shown by the blue dashed 
line. Nine out of 33 objects would fall within our i-dropout selec- 
tion. Hence, the i-band dropout technique seems to select Ri30% 
of the z s=s 5.7 Lya emitters found by narrowband surveys. It is 
relevant to note that four of these ar e the weakest emitters of the 
confirmed sample of Hu et al. (2001). The other objects do not 
meet our selection criterion because of the contribution of Lya flux 
to the i-band magnitude (see Figure^). As noted in Fig. [?] this 
selection bias is not a problem for higher redshift (z > 6) LAEs. 



31% ar e found by Stanwav et al. ( 2004a) , Stanwa v et al.l 
(2004bJ), and I Vanzella et all (|2006h . These observations, 
and the results of the Lya emission completeness test of 
£|4.3i suggest that Lya emitting objects represent «30% 
of i-dropout spectroscopic samples. 

However, just because the fraction of i-dropouts with 
Lya emission is 30% does not imply that the same thing 
is true for the 2^6 population as a whole. In tj4.31 
and easily seen in Fig. [JJ the 4775 — 2350 >1.3 color cut 
biases our selection against the inclusion of lower redshift 
(z < 5.8) objects with strong Lya emission (>50A). This 
suggests that the effective selection volumes (and redshift 
distributions) of i-dropouts with significant Lya emission 
may be quite different from the distribution without such 
emission. 

This can be acc o unted for by utilizing the simulations 
of IBouwens et all (|2006l ) (see Fig. [7J to determine the 
effective search volumes, and corresponding densities, 
for objects of different Lya equivalent widths. For the 
distribution of Wh ya in our sample (i.e., three objects 
with Wlj/q=25 A, two with H / Ly a =50 A, and one with 
WLya— W0 A) these simulations imply that 32 ± 10% of 
star-forming galaxies at z ~ 6 have WL ya >20 A, and 
7 ± 6% have W Lya >100 A. 12 In the situation where 

12 We compute the errors on the fraction of LAEs assuming a 
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every Lya emitting object in our sample has Wl v 



-100 



A and spectroscopy misses every object with z > 6.f5, 
this fraction increases to 46 ± 11%. This case is an up- 
per limit given our assumptions. A lower limit of 25 ± 
10% is determined, if the emission line objects GOODS 
i6 and UDF PFs i4 are lower redshift contaminants. A 
similar Lya emitting fraction of 32 ± 14 % (upper limit 
of 47 ± 16%) is found using the results of lStanwav et al.l 
(|2004afi . 

These results are similar to those obtained from 
Shapl ev" 3t al.1 (|2003l ), who analyzed spectr a of over 1000 
z ~ 3 LBGs. Though [Shaplev et alj (|2003[ ) selected their 
LBGs photometrically using a two-color cut, the extra 
color is used to remove potential low-redshift contami- 
nants. As noted in § 14.11 the expected contamination 
of our survey is small; hence, we believe our z ~ 6 
LBG survey i s simil ar enough to warrant comparison. 
IShaplev et all (|2003ft find that only 25% of their sam- 
ple of z ~3 LBGs have W Lya >20 A, and -2% have 
WL ya >100 A. This is similar to the value we deter- 
mine from our observe d Lya fraction (as well as that of 
IStanwav et alj (|2004a| )). Moreover, the fraction of z ~ 
3 objects with Wi, ya >20 A is within 2a of our z ~ 6 
upper limit. Therefore, it appears that there is no strong 
evolution in the fraction of Lya emitting objects between 
z ~ 3 to z ~ 6. 

Differing results are found by Shimasaku ct al. ( 2006), 
who compare the z ~ 6 LBG luminousity function of 
iBouwens et al.1 (|2006l ) with that derived from their z w 
5.7 LAE sample. These authors determine that nearly 
every z ~ 6 LBG with M uv < - 20 (z 8 50 < 26.6 
mag) should have Wl V o >20 A. In fact. lShimasaku et al.1 
(2006) argue that the fraction of z ~ 6 LBGs with Wr, ya 
>100 A is «80%. Clearly these result s are inconsistent 
with o ur analysis above. At best, the [Shimasak u et al.1 
(2006) fraction is a factor of «2 too high. To be con- 
sistent with our survey, the spectroscopic efficiency we 
calculated in § l4.3l must be a factor of >10 too low, a sit- 
uation we consider highly un l ikely given our agreement 
with th e results of iHu et all (|2004h and IStanwav et al.1 
(|2004af) . Hence, we believe tha t the Lya emitting frac - 
tion of nearly unity found by fShimasa ku et al.1 (2006) 
must be wrong. This may be at least partially due 
to the mild evolution in the luminosity function which 
will occur between z w 5.9 (the mean redshift of the 
IBouwens et al.1 (|2006l ) i-drop out selection) to z ~ 5 .7, 
which is not accounted for by IShimasaku et alj ((2006) . 

5.2. Luminosities and Star Formation Rates 

Objects in our sample with measurable Lya emission 
all have observed Z85o-band magnitudes of 26.1 or fainter, 
populating the faintest 70% of the upper panel of Fig- 
ure [6] Furthermore, the brightest two objects in our 
CL1252 sample were confirmed z ~ 6 LBGs without any 
detectable Lya emission. However, after delensing the 
sources in our CL1252 sample (lower plot of Fig. O, 
this luminosity segregation of Lya emitti ng objects is re- 
moved. The Lya emitting LBGs found by Sta nwav et al.l 
(|2004al) have z 850 -band magnitudes of 25.48 ± 0.03, 26.37 
± 0.06 , and 27.15 ± . 12. Asi de from the brightest object 
in the IStanwav et a l. (2004a) sample and the extremely 

binomial distribution. This assumes systematic errors are small. 



bright (.2850 = 24.7) object of lBunker et~aTl (|2003|L all 
confirmed z ~ 6 Lya emitting z-dropout objects have 
zsso-band magnitudes in excess of 26. 

The average star formation rate derived from the UV 
continuum (SFRuy) of our entire sample (excluding the 
6L* object BD38) is 8.8 ± 3.9 h^ 2 7 M vr -1 . This value 
increases slightly to 9.2 ± 4.4 M Q yr^ 1 if we exclude Lya 
emitting objects. The average SFRjjv of Lya emitting 
objects is 7.7 ± 2.2 M© yr _1 . Hence, the SFR UV of 
Lya emitting objects appears to b e similar to tha t of the 
i-dropout population in general. IShaplev et al.1 (2003) 
find at z ~ 3 that Lya emitting galaxies with stronger 
Lya emission have lower star formation rates. However, 
as noted by these authors, this trend could be at least 
partially due to selection effects. Given the relatively 
small size of our sample, and the large uncertainties (up 
to a factor of 10) inherent in UV continuum derived SFRs 
(see iPapovich et~aT1l2005h , our results do not allow us to 
make strong statements about the SFRs of Lya emitting 
LBGs relative to the population as a whole. 




SFR 



UV 



Fig. 9. — Histograms of SFRjjv for our spectroscopic sample. 
Sources in the CL1252 field were corrected for lensing by the clus- 
ter potential, except for objects BD22 and BD62 where the cor- 
rection was highly uncertain. The solid filled histogram are Lya 
emitting objects, and the cross hatched are the confirmed z ~ 6 con- 
tinuum objects lacking detectable Lya emission. The Lya emitting 
objects appear to have SFRuy very similar to the sample average 
(excluding the 6L* galaxy BD38). 



5.3. Correlation between Lya Emission and Galaxy Size 

In Figure [TU] we show zsso-band images (l'/5 thumb- 
nails) for our entire spectroscopic sample (excluding 
UDF PFs iO). The uppermost six thumbnails are the 
confirmed Lya emitting LBGs, the two thumbnails in 
the bottom row are the confirmed z ~ 6 objects with 
no detectable Lya emission, and the middle 13 are ob- 
jects with no clear spectral features. The Lya emitting 
LBGs all have compact morphologies and show no mor- 
phological disturbances. However, the morphologies of 
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Fig. 10. — Thumbnail images of the entire spectroscopic sample in the zgso-band (l'J 5 on a side). UDF PFs is abbreviated as UPF. 
The upper six thumbnails are our Lyct emitting sample, the lowest two are the confirmed z ~ 6 objects that are undetected in Lya (BD38 
and BD03), and the middle 13 are objects without any observed spectral features. The morphologies of BD38 and BD03 are much more 
disturbed than the Lya emitting LBGs. Furthermore, the majority of the undetected, fainter objects have a similar distorted morphology, 
in contrast to the very compact appearance of the Lya emitting LBGs. 



BD38 and BD03 are very disturbed and extended, with 
"plumes" extending from a "core" . Furthermore, the un- 
detected objects overall have much more extended mor- 
phologies tha n the objects w i th dete cted Lya emission. 

Similarly, iStanwav et al.l (|2004al lbl) confirmed Lya 
emitting z-dropout objects have compact morphologies 
in the z 850 band, with 0'.'09 < r w < 0'/14. More- 
over, SBM03#01 also appears compact in HST NIC- 
MO S 1.1 Mm FllOW-b and and 1.6/x F160W-band imag- 
ing (|Evles et al.ll2005[) . Interestingly, the Lya-emitting 



z = 5 78 galaxy SBM03#03 discovered bv lBunker et~aTl 
(2003) also has a very compact morphology (r^ < 0'.'08; 
marginally resolved in ACS imaging), despite being ex- 
ceedingly bright (z 850 = 24.7) for an i-dropout. This is 
in contrast to BD38, a similarly bright (zsso = 24.6) i- 
dropout in our spectroscopic sample, but which is much 
more extended (tm — 0'.'25) and lacking Lya emission. 

In the upper panel of Figure [Tl] we plot the thi 
versus zsso-band magnitude for sources in our spec- 
troscopic sample. The red circles are Lya emitting 
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Fig. 11. — (Upper -panel) Half-light radius versus zgso-band mag- 
nitude. The small purple dots rep resent a sample of 332 i-dropouts 
selected over the GOODS fields (Bouwens et al\ \200b). The red 
stars a re previously spectroscopically- confirmed z ~ 6 Lya emitting 
LBGs (Bunker et al. 2005; Stanwav et al. 2004a b; Vanzella et al. 
2006,). The other symbols are the same as FiaureW\ The black dot- 
dash-dot-dot li n e is t he surface brightness selection limit of the 
I Bouwens et al\ $2006) sample. The purple line denotes the best 
linear fit to this sample (purple dots; other objects excluded). The 
vertical orange error bar of a = 0.05 dex is the scatter in the fit. 
The horizontal error bar is the typical photometric uncertainty. 
The blue dash-dot line is the FWHM of the zgso-band PSF and 
represents a reasonable floor below which our size measurements 
become g uite uncertain . (Low er panel) Residuals from the linear 
fit. The I Bunker et ali fB005\ ) object is outside the plot window 
(Alogrf,; = — 0.49). All Lya emitting objects hav e half-light radii 
smalle r than is typical for the 332 i-dropouts of Bouwens et al. 
(200b ). A Rank-Sum test performed on the objects from our spec- 
troscopic sample supports the hypothesis (at >97% confidence) that 
Lya emitting i-dropouts (6 objects from our sample) are smaller 
than the objects without such emission ( 14 objects from our spec- 
troscopic sample). Furthermore, the Lya emitting objects have a 
fairly constant size (r^i = 0.11 ± 0.01) over the entire range of 
luminosities probed. 



LBGs, the blue triangles represent BD38 and BD03 
(sources with a continuum break), and the black squares 
are i-dropouts in our sample without any clear spec- 
tral features. The CL1252 objects have been delensed. 
Sh own with purple dots a re the 332 i-dropouts identified 
by iBouwens etafl (|2006f ) over the two GOODS fields. 
Also included (red stars) are the four prev i ously con- 
firmed Lya em i tting L BGs of iBunker et al.l (|2003f ) and 
iStanwav et alj (|2004al lbl) noted above, where the half- 
light radii and Zssn-band magnitudes quoted are taken 
from IBouwens et al.l (|2006Tj 13 . The purple line in Fig. [Til 
de notes a linear fit to the i-dropout sample (purple dots) 
of IBouwens et all ()2006l ). Objects with half-light radii 
smaller than the FWHM of the z 850 -band PSF (0"09; 
blue dash-dot line in figure), fainter than our magnitude 
limit (zs50 > 27.8), or brighter than z 850 = 25.0 (due to 
a scarcity of bright sources) were not included in the fit. 



13 Object GOODS-N i'Drop 6 of IStanwav et alj |2004a) is not 
included, due to its uncertain line identification 



No sources from our spectroscopic sample were included 
in the fit. 

We want to test the hypothesis that sources with Lya 
emission are smaller on average than the i-dropout pop- 
ulation in general. We chose not to use a \ 2 test, because 
it assumes a Gaussian scatter around the best-fit distri- 
bution. A better test is the Rank-Sum test utilizing only 
the half-light radii of our spectroscopic sample. Using 
this test to compare our spectroscopic sample with the 
IBouwens et alj (|2006f) i-dropout sample, we find that the 
half-light radii of the galaxies in our spectroscopic sam- 
ple are s ignificantly (at > 94% confidence) smaller than 
those of IBouwens et all ( 20061 ). Restricting the Rank- 
Sum test to include only objects in our spectroscopic 
sample also supports the hypothesis that Lya emitting 
i-dropouts are smaller than the rest of our sample at 
>97% confidence. Hence, Lya emitting i-dropouts seem 
to be morphologically distinct from other i-dropouts. 

One p otential except i on to this rule is the object UDF 
5225 of iRhoads et~atl (l2005h . which has a "plume" of 
1" extending from a compact "core". This object was 
selected as a Vesoe-dropout, has a spectroscopic redshift 
of z = 5.480, and is too blue (1775 — z 850 « 0.5) to be 
selected by i-dropout methods. However, even if we as- 
sumed that this object passed our color selection, the 
surface brightness of the plume is below our limiting sur- 
face brightness, due to shallower photometry than the 
UDF. Our deepest photometry (UDF PFs) of the object 
would only contain the brighter compact "core" . Hence, 
this object would not appear abnormal for Lya emitting 
objects in our sample. 

Note that the nature of the vast majority of the ob- 
jects (332 in total) used in our larger i-dropout sample 
is not known; however from the discussion above, ap- 
proximately 30% of these objects should be Lya emitting 
LBGs. If so, this would mitigate the differences observed 
between the sizes of Lya emitting LBGs and those from 
the i-dropout population in general (which include a sig- 
nificant fraction of Lya emitting sources). Hence, the 
size-luminosity discrepancy between Lya emitting and 
non-emitting LBGs would be larger than quoted here. 

This morphological difference may be due to the in- 
fluence of the Lya emission on the zsso-band flux. To 
date, the only confirmed z ~ 6 Lya emitting i-dropou t 
with NICMOS imaging is SBM03#01 (lEvles et al.ll2005h . 
which, as noted above, is compact in both 1.1 /im 
FllOW-band and 1.6/mi F160W-band imaging. Hence, 
the compact morphology of this LBG cannot simply be 
the result of its Lya emission alone. NICMOS imaging 
of many more i-dropouts with confirmed Lya emission 
is needed to determine the overall relevance of Lya emis- 
sion to the Z850-band morphology of these objects. 

One potential physical explanation of this size devia- 
tion could lie in the masses of these objects. Lya emitting 
objects at z ~ 3-4 are found to have smaller stellar masses 
than objects lacking this emission (jOverzier et alj|2006t 
iGawiser et alJl2006t ). This suggests that these objects 
may be associated with less massive dark matter haloes. 
One would make a similar inference for the masses of 
these objects using th e observation that these objects are 
less d usty on average (IShaplev et al.ll2003t IGawiser et al.l 
2006), thereby having fainter dust-corrected luminosities, 
than sources with no observable Lya emission. Assum- 
ing this translates to z ~ 6, and using the well-known 
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correlation between size and mass (|Mo et al.lll998h . we 
would expect Lya emitting objects to appear smaller on 

average, which is what we observe. 

We note that the recent results of lLai et~a l. (2006) on a 
sample of 12 z ~ 5.7 LAEs over the HDF-Nort h GOODS 
area c ould be seen to support this finding. lLai et all 
(2006) found that 3 of the 12 sources were clearly de- 
tected in the IRAC data and inferred masses of 10 9 Mq 
to 10 10 Mq, through detailed stellar population mod- 
elling. While these masses are comparable to similar 
luminosity i-dropouts studied ove r the GOODS fields 
(|Yan et al.lf2006l: lEvles i et alll2006l ). the vast majority of 
LAEs in the lLai et al.l ( 20061 ) sample (9 out of 12) are not 
detected in the GOODS IRAC imaging and hence will 
be significantly less massive than the 3 IRAC-detected 
LAEs. Though clearly a more careful comparison be- 
tween these populations is needed, this suggests that 
z ~ 5.7 LAEs, on average, are somewhat less massive 
than the typical z ~ 6 star-forming galaxy. 

5.4. Preselecting Galaxies with Lya Emission 

One interesting aspect of this morphology dependence 
on Lya emission strength is in its potential to prese- 
lect Lya emitting i-dropouts for spectroscopic follow-up, 
similar to what is achieved in narrowband surveys. Con- 
firmation rates of i-dropout selected z ~ 6 LBGs are 
rather low («30%) due to the low S/N of these objects 
in the continuum and the small fraction of LBGs with 
Lya emission (§5.1). From our spectroscopic sample, we 
note that a simple size cut Alog(rhi) < -0.05 and magni- 
tude cut z 850 > 25.9 would select five (83%) of the Lya 
emitting LBGs, while rejecting all but two (15%) of the 
sources without detectable Lya emission. This would in- 
crease our success rate for spectroscopic confirmation by 
more than a factor of three. 

6. CONCLUSIONS 

We have obtained spectroscopic observations of 22 
i-dropouts drawn from two deep ACS fields. These 
dropouts were selected to have 1775 — zsso colors greater 
than 1.3. Spectroscopic redshifts for eight z ~ 6 objects 
in the RDCS 1252-2927 and Ultra Deep Field Parallel 
fields were derived. We thereby confirm the eff e ctive- 
ness o f the i-dropout technique of iBouwens et all (2003, 
l2004a| ) in isolating a statistically relevant sample of 5.5 
< z < 6.2 star-forming galaxies. 

No clear case of contamination by low-redshift sources 
was found in this sample. Together wit h the results of 
the co mplementary ACS Grism survey oflMalhotr a et all 
( 2005), this suggests our spectroscopic sample is domi- 
nated (>90%) by galaxies at z ~ 6. Six of the confirmed z 
~ 6 objects possess measurable Lya emission, with z S50 - 
band magnitudes ranging from 26.1 to 27.1. The two 
brightest objects in our CL1252 sample are continuum 
sources, lacking detectable Lya emission, but show clear 
evidence for a Lyman break. 

We compare our findings with other z ~ 6 i-dropout 
surveys, z w 5.7 narrowband surveys, and with the global 
properties of LBGs at z ~ 3 to determine the nature 
and composition of z ~ 6 galaxies. Our findings are as 
follows: 

1. Significant (W Lya > 20 A) Lya emission is de- 
tected in the spectra of only 30% of i-dropout ob- 
jects. Utilizing the model redshift distributions of 



IBouwens et al.l (|2006| ) to control for selection bi- 
ases, we infer that only 32±10% of star-forming 
galaxies at z ~ 6 show significant (WL ya > 20A) 
Lya emission (§5.1: Figure [7]), with an upper 
bound of 46±11%. Moreover, the Wh ya of these 
objects on average are much smaller than those 
found for narrowband-selected LAEs. Since these 
tren ds are also evident in the LBG population at z 
~ 3 (jShaplev et al.ll2003l ). this suggests that there 
is no strong change in the fraction of Lya emitting 
objects in LBG population from z > 5 to 2 < 4. 

2. The i-dropout technique misses ~70% of narrow- 
band selected z rs 5.7 LAEs. Moreover, the LAEs 
i-dropout surveys do select show the weakest Lya 
emission. This is a consequence of the Lya emis- 
sion in the i775-band flux at z < 5.9. However, 
the selection efficiency increases strongly with red- 
shift, due to Lya emission shifting from the 2775 to 
Z850 bands. By z ~ 6, every zgso-band detected 
Lya emitter would be sampled by z-dropout tech- 
niques. Using the simulations in £15.11 we expect 
to miss 10-46% of the LAE population with our 
i-dropout selection. 

3. Lya emitting LBGs have similar zgso-band mag- 
nitudes and UV SFRs to those without detectable 
emission. However, no Lya emitting object in our 
sample has zsso < 26.1. 

4. The size of Lya emitting objects is more compact 
than predicted from the observed size-luminosity 
relation of 332 i-dropout galaxies. One possible 
explanation for this trend is that sources with Lya 
emission may be systematically less massive than 
the typical i-dropout. 

Most importantly, these results suggest that Lya emit- 
ting objects only constitute a modest fraction of the 
LBGs at z ~ 6 and thus do not provide us with a repre- 
sentative sample. 

These results have important implications for how sur- 
veys for star-forming galaxies at z ~ 6 should be con- 
ducted. First, it appears that our spectroscopic suc- 
cess rates can be significantly enhanced by targetting 
i-dropouts with compact zsso-band morphologies. From 
our sample, we find the criteria zsso > 25.9 and A logr/,.; 
< -0.05 selects >80% of Lya emitting objects, while re- 
jecting >80% of the remaining sample. 

Second, it appears wide-area narrowband surveys miss 
«70% of the broadband-selected LBG population at z ~ 
6, and as seen in FigurefTTIthe objects they do find are in- 
trinsically different from the i-dropout LBG population 
as a whole. On the other hand, broadband surveys miss 
the strongest Lya emitters at 5.5 < z < 5.8 that nar- 
rowband surveys sample. Therefore, it appears a combi- 
nation of surveys - taking advantage of the strengths of 
both narrowband and broadband selection techniques - 
will be necessary to obtain a complete characterization 
of the star-forming population at z ~ 6. 
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Fig. 2. — LRIS spectra and ACS imaging of emission line objects: two-dimensional unsmoothed spectra (upper panel), candidate Lyct 
emission (lower panel), and ACS imaging (lower postage stamps). Also shown is the two-dimensional a=l pixel Gaussian smoothed spectra 
of the sky (middle panel). Objects are (a) BD46, (b) BDOO, (c) GOODS 16 0, (d) UDF PFs H, (e) UDF PFs it, and (f) UDF PFs %B. The 
postage stamps are l'.'O X l'.'O. All upper limits given for the B453, Veo6j an d *775 band fluxes are 2 a. The exposure times for the objects 
were 7200 s, except for the two CL1252 objects (BD46 and BDOO) where the exposure times were 16200 s. Each one- dimensional spectrum 
has been smoothed with a 3 pixel boxcar filter, except BDOO and BD46. Vertical dotted lines delineate the region used for measuring line 
fluxes, equivalent widths, and FWHMs. For the UDF PFs objects, the blue cross hatched rectangles represent regions where strong skylines 
are present. The 1775 — 2850 colors and the asymmetry of these line profiles are consistent with the emission lines from BD46, BDOO, UDF 
PFs il, and UDF PFs i2 being Lyct X1215.67. The other two emission lines (GOODS i6 and UDF PFs H) are difficult to identify, due 
to skyline interference. We consider the implication that the features in these latter two sources are spurious and Lyct is undetected. 
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Fig. 3. — (a) FORS2 two- and one dimensional spectra of BD46 (LRIS spectrum shown in Fig. 0). The total exposure time is 22.3 
hours. The Lyce emission occurs at 84 06A, consistent with the LRIS spectrum, (b) Shows the faint continuum of BD46 redward of the 
emission. The two dimensional spectrum has been smoothed with a cr=l gaussian. The one dimensional spectrum has been smoothed with 
a 10-pixel boxcar filter. The red dashed line delineates the continuum of f = 3.4 ± 1.0 counts A^ 1 . 
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Fig. 4. — FORS2 spectra and ACS/ISAAC imaging of BD03 (z = 5.98 ± 0.1). The upper two dimensional spectrum is unsmoothed, 
while the lower one is smoothed with a a=l pixel Gaussian. The extracted one- dimensional spectrum is smoothed with a 5 pixel boxcar 
filter. From left to right, the lower four images are irrs, 2850; J, o/nd K. The total exposure time is 22.3 hours. The lower axis is the 
observed wavelength shifted to the systematic redshift z = 5.98. Blue arrows indicate the position of the continuum break. The dotted lines 
represent the wavelength region used to fit the continuum near the break, where error bars represent the 1 a error. 
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Fig. 5. — FORS2 spectra, ACS and ISAAC imaging of BD38: two-dimensional spectrum and extracted one- dimensional spectrum 
smoothed with a 5 pixel boxcar filter. From left to right, the lower four images are 1775, 2:850, J, and K. The total exposure time is 22.3 
hours. The lower axis is the observed wavelength shifted to the systematic redshift z = 5.515. Blue arrows indicate the position of the 
continuum break. The dotted lines represent the wavelength region used to fi t the continuum near the break, where error bars represent 1 
a error. A more complete spectrum is shown in Dow-Hyqelund et al. (2005). 



